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> ; ABSTRACT 

. We have constructed a 3D photoionisation model of a planetary nebula (PN) similar 

' in structure to NGC 7009 with its outer pair of knots (also known as FLIERs -fast, 

^ , low-ionization emission regions). The work is motivated by the fact that the strong [N 

I ■ li]A6583 line emission from FLIERs in many planetary nebulae has been attributed 

ly-^ ' to a significant local overabundance of nitrogen. We explore the possibility that the 

, apparent enhanced nitrogen abundance previously reported in the FLIERs may be due 

' to ionization effects. The model is constrained by the results obtained by Gongalves, 

^ ; Corradi, Mampaso & Perinotto (2003) from the analysis of both HST [O ill] and [N 

Y^, li] images, and long-slit spectra of NGC 7009. Our model is indeed able to reproduce 

Q ' the main spectroscopic and imaging characteristics of NGC 7009 's bright inner rim 

^ I and its outer pairs of knots, assuming homogeneous elemental abundances throughout 

C/3 ■ the nebula, for nitrogen as well as all the other elements included in the model. 

' We also study the effects of a narrow slit on our non-spherically symmetric density 

^ , distribution, via the convolution of the model results with the profile of the long-slit 

• '~j ■ used to obtain the spectroscopic observations that constrained our model. This effect 

' significantly enhances the [N ii]/H/3 emission, more in the FLIERs than in the inner 

rim. 

Because of the fact that the (N+/N)/(0+/0) ratio predicted by our models 
are 0.60 for the rim and 0.72 for the knots, so clearly in disagreement with the 
N"'"/N=0"'"/0 assumption of the ionization correction factors method (icf), the icjs 
will be underestimated by the empirical scheme, in both components, rim and knots, 
but more so in the knots. This effect is partly responsible for the apparent inhomo- 
geneous N abundance empirically derived. The differences in the above ratio in these 
two components of the nebula may be due to a number of effects including charge 
exchange -as pointed out previously by other authors- and the difference in the ion- 
ization potentials of the relevant species -which makes this ratio extremely sensitive 
to the shape of the local radiation field. Because of the latter, a realistic density dis- 
tribution is essential to the modelling of a non-spherical object, if useful information 
is to be extracted from spatially resolved observations, as in the case of NGC 7009. 

Key words: Atomic data - ISM: abundances - planetary nebulae: individual 
(NGC 7009). 



1 INTRODUCTION 

It is well known that planetary nebulae (PNe) possess a 
number of small-scale structures that, at variance with their 
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large-scale structures, are prominent in low-ionization emis- 
sion lines, such as [N ll] and [O ll]. In terms of morphology 
and kine matics, the low-io nization structures (LISs, for a re- 
view see iGoncalvejboO-d) may appear as knots, filaments, 
jets, and isolated features that, in some cases, move with 
the same velocity of the ambient nebula, or supersonically 
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Table 1. Main pliysical/cliemical parameters from Paper I. 



Parameter R K NEB 

Eastern Western Eastern Western Whole Nebula 



Ne[S ll](cm-3) 


5,500 


5,900 


2,000 


1,300 


4,000 


Ne[Cl m](cm-3) 


5,200 


5,900 




1,900 


1,300 


Te[0 m](K) 


10,000 


10,200 


9,600 


10,400 


10,100 


Te[N n](K) 


10,400 


12,800 


11,000 


11,700 


10,300 


Te[S n](K) 






7,100 


9,400 




He/H 


1.08{-1) 


1.16(-1) 


1.02(-1) 


9.55(-2) 


l.ll(-l) 


O/H 


4.5(-4) 


4.82(-4) 


5.8(-4) 


4.5(-4) 


4.71(-4) 


N/H 


7.0{-5) 


1.8(-4) 


3.8(-4) 


2.5(-4) 


1.7(-4) 


Ne/H 


l.l{-4) 


l.l(-4) 


l.l{-4) 


1.3(-4) 


l.l(-4) 


S/H 


6.1{-6) 


4.9(-6) 


1.39{-5) 


9.3(-6) 


8.3(-6) 



Empirically derived parameters, for the different structures of the NGC 7009, for both, the Eastern and Western sides of the nebula, 
with "R", "K" and "NEB" standing for rim, outer knots and the whole nebula (see Tables 1 and 3 of Paper i). 



through the envi ronment, like FLIERs (fast, low- ionization 
emission regions. lBalic^^^lJ^ 93 | 'l or BRETs (bipolar, ro- 
tating, episodic jets : lL6pez et alJll995^ . 

For over ten years, the strong [N ll] line emission 
from LISs has been attributed to a significant local over- 
abundance of nitro gen (Balick et al. 19 93, 1994, 1998; and 
references therein'!. iBalick et al.l (|^9^ interpreted the N- 
enrichment in FLIERs as evidence of their origins in recent 
high- velocity ejections of the FN central star. The above 
work and successively IBalick et alJ lll994ll , which included 
the derivation of the ionic and elemental abundances of 
LISs, pointed out "an apparent enhancement of nitrogen 
relative to hydrogen by factors of 2-5" in FLIERs. How- 
ever, since th en, further work casted doubts over the lat- 
ter statement ian et al .ll997l:lAlex ander fc Bahckll997t 
IGoncalves et alJl2003t IPerinotto et al . 2004a). 

Abundances in PNe can be derived from the analysis 
of coUisionally excited line (CEL) or optical recombination 
line (ORL) spectra using empirical metho ds or tailored pho- 
toionisation models (e.g. IStasil^skall2004^ or a combination 
of the two. For some elements (e.g. He) and ions (e.g. 
if only the optical spectrum is available) the use of ORL is 
mandatory, whilst for others the use of CELs may be un- 
avoidable. But, when both CEL and ORL determinations 
are available, the former have been commonly preferred be- 
cause they are generally stronger and easier to detect than 
ORLs. Moreover, there is a well known discrepancy between 
CEL- and ORL-abundanccs, as well as electron temperature 
determinations (Liu ot al. 1995). Empirical abundance anal- 
ysis rely on ionization correction factors (icjs) to account for 
the unseen ions (e.g. Kingsburgh & Barlow, 1994). Results 
obtained with the icf method can be somewhat uncertain 
in some cases, particularly when they arc applied to spa - 
tially resolved long-slit spectra (Alexa nder fc Balickl | l997| ). 
as has been t he ca se for the work of ba lick et aLr i^94li ; 
iHaiian et alJ il997tl and IGoncalves et alJ ~ (l200^ ~hereafter 
Paper i). 

In fact, the abundances derived in Paper l and shown 
here in Table 1, from optical long-slit spectra of NGC 7009 , 
using the icf scheme of IKingsburgh fc Barlowl il994l) . 

showed only a marginal evidence for overabundance of N/H 



in the outer knots of the neb ula (the ansae), rein forcing the 
doubts over previous results jBalick et al.lll994^ where the 
N/H enhancement of a factor of 2-5 in the ansae were re- 
ported. 

One of the major shortcomings of empirically- 
determined chemical abundances lies in the fact that a num- 
ber of assumptions on the ionization structure of the gas 
need to be made in order to obtain the icfs. A preferred 
alternative would be the construction of a tailored photoin- 
isation model for a given object, aiming to fit the emission 
line spectrum and, in the case of spatially resolved objects, 
projected maps in a number of emission lines. 

NGC 7009, the "Saturn Nebula", is a PN compris- 
ing a bright elliptical rim. Its small-scale structures in- 
clude a pair of jets and two pairs of low-ionization knots. 
On a larger scale, it is known that NGC 7009 possess 
a tenuous hal o with a diameter of more than 4 arcmin 
iMoreno et alj |199 8), whose inner regio ns display a system 
of concentric rings jCorradi et alJl2004l) like those observed 
in NGC 6543 and other few PNe llBalick et ajbOOlfl . High- 
excitation lines dominate the inner regions along the mi- 
nor axis, while emission from low-ionization species is en- 
hanced at the extremities of the major axis. The ioniza- 
tion structure is further enriched by the fact that the low- 
excitation regions present strong variations in excitation 
level and dumpiness. NGC 7009 was classified as an oxygen- 
rich PN llHvung fc AlleJll995ah . with an O/C ratio exceed- 
ing 1, and anomalous N, O, and C abund ances iBakeilll983l : 
IBalick et al]ll994l : iHvung fc AlleJll995a^ . Its central star is 
an H-rich 0-type star, with effective temperature of 82 000 K 
( Mcndcz ct al. 1992; Kingsburgh & Barlow 1992) . The kine- 
matic s of NGC 7009 was stud ied first bv lReav fc AthertonI 
(^8^ and IBalick et al.l il987fl . who showed that the ansae 
are expanding near the plane of the sky at highly super- 
sonic velocities. The derived inclination of the inner (caps) 
and outer (ansae) knots, with resp ect to the line of sight, ar e 
i = 51° and i = 84°, respectively llReav fc Athertonlll985ll . 
More recently, IFernandez et alj i2004ft have measured the 
proper motion and kinematics of the ansae in NGC 7009, 
assuming that they are equal and opposite from the central 
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star, obtaining Vcxp ~ 114 ± 32 kms ^, for the distance of 
- 0.86 ± 0.34 kpc. 

In this paper we will present a simple 3D photoionisa- 
tion model, aiming at reproducing the observed geometry 
and spectroscopic "peculiarities" of a PN like NGC 7009, 
exploring the possibility that the enhanced [N ll] emission 
observed in the outer knots may be due to ionization ef- 
fects. We wil l use the 3D photoinisation code MOCASSIN 
jErcolano et al. 2003a'l together with the long-slit data pre- 
sented in Paper I and the HST images available for this 
nebula. Our model and the observational data are described 
in Section 2. The results for our main model are presented 
and discussed in Section 3, while in Section 4 we discuss an 
alternative model, aimed at highlighting the relevance of the 
geometry and density distribution in such a kind of mod- 
elling. The work is briefly summarised in Section 5, where 
our final conclusions are also stated. 



2 METHOD 

2.1 Observational data 

The observations used to constrain the photoionisation 
model were described in detail in Paper I. These included 
HST [O III] and [N ll] images as well as Isaac Newton Tele- 
scope long-slit, intermediate dispersion spectra, along the 
PN major axis (P.A. = 79°). For further details see Paper I. 

We refer to Figure 1 of Paper I, the HST [O ill] and 
[N ll] images of NGC 7009, in which the nebular features 
are outlined. Here we will employ the same nomenclature 
introduced in Paper I. See also the left panel of Figure 2, in 
Section 3.5. 

2.2 The 3D photoionisation code: MOCASSIN 

The nebula w as modeled using the 3 D photoinisation code, 
MOCASSIN, of lErcolano et al.l Il2003aih . The code employs a 
Monte Carlo technique to the solution of the radiative trans- 
fer of the stellar and diffuse field, allowing a completely 
geometry-independent treatment of the problem without 
the need of imposing symmetries or approximations for the 
transfer of the diffuse component. 

The reliability of the code was demonstrated via a 
set of benchmarks describ ed by iPeouignot et alJ (1200 ll) 
and lErcolano et al.l (|20033). A number of axy-symmetric 
planetary nebulae have already be en modeled using MO - 
CASSIN, exa mples include NGC 39 18 iErcolano et all2003bl) . 
NGC 1501 tercolano et al."2004') and the H-deficient knots 
of Abell 30 fercolano ct al. 2003al . 

2.3 Input parameters 

A thorough investigation of the vast parameter space was 
carried out in this work. This involved experimenting with 
various gas density distributions, central star parameters 
and nebular elemental abundances. The model input pa- 
rameters that best fitted all observational constraints are 
summarised in Table 12.31 and discussed in more detail in 
the following subsections. 

A common problem when studying galactic PNe is the 
large uncertainties associated with the distance estimates. 



Table 2. Input parameters for the model. 



L* {Lq) 


3136 


N/H 


2.0(-4) 


'^eff (K) 


80,000 


O/H 


4.5(-4) 


Rin (cm) 


0.0 


Ne/H 


1.06 (-4) 


Rout (cm) 


3.88(17) 


S/H 


0.9(-5) 


He/H 


0.112 


Ar/H 


1.2(-6) 


C/H 


3.2(-4) 


Fe/H 


5.0(-7) 



Abundances are given by number, relative to H. 

that propagate to the determination of the nebula geom- 
etry and central star parameters. We adopted a distance 
of 0.86 kpc for our mo dels of NGC 7009, as computed by 
iFernandez et al .' ( 2004) from the weighted average of 14 val- 
ues iAcke^^l. 1992) determined with statistical methods; 
the value was quoted with an error of ± 0.34 kpc. 

2. 3. 1 Stellar parameters 

After having experimented with various stellar atmosphere 
models to describe the ionising continuum, we reverted to 
using a blackbody of T^// = 80,000 K and log L* = 3.50, 
as this res ulted in the best fit of the nebular emission line 
spectrum. iMendez et al.l lll992l) determined the Te// and 
log g for the central star of NGC 7009 using non-LTE model 
atmosphere analysis of the stellar H and He absorption line 
profil es, finding Tgf f = 82 000 K and log L« = 3.97, in solar 
units. iMendez et all il992l) assumed a distance of 2.1 kpc in 
their analysis which therefore resulted in their value for the 
stellar luminosity being higher than that inferred from our 
modelling. 

2.3.2 Elemental Abundances 

The nebular elemental abundances used for the photoioni- 
sation model are listed in Table 2, where they are given by 
number with respect to H. Although the MOCASSIN code can 
handle chemical inhomogeneities, these were not included in 
our models as they proved to be not necessary to reproduce 
the GEL spectra of the R and K regions of NGC 7009. 

The values shown result from an iterative process, where 
the initial guesses at the elemental abundances of He, N, O, 
Ne and S, ta ken from Paper I (see Table 1), and those of C 
and Ar, from lPottaschI ll2000l) . were successively modified to 
fit the spectroscopic observations. 

2.3.3 Density distribution 

The simplest possible density distribution model was con- 
structed in order to demonstrate that the spectroscopic pe- 
culiarities often found in LISs can be the product of sim- 
ple (and well-known) photoionisation effects. The case of 
NGC 7009 is taken as an example, but the emphasis is not 
in the construction of a detailed model for this object in 
particular. With this in mind we described the nebula by 
an ellipsoidal rim with a H number density, Nh, peaking 
to 9000 cm"'' in the short axis direction exponentially de- 
creasing to a minimum value of 4000 cm~^ in the long axis 
direction. The short and long axes of the inner and outer 
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ellipsoids measure 3.84x10^ cm and 9.99x10^ cm, and 
7.06x10^® cm and 1.84x10^'^ cm, respectively, at the dis- 
tance assumed for NGC 7009. The rim is surrounded by a 
spherical shell of less opaque, homogeneous density gas, with 
Nh = 1600 cm~'^. The diameter of the sphere is equal to the 
long axis of the outer ellipsoid defining the rim. Cylindrical 
jets, 1.75x10'^'^ cm in diameter, connect the rim to a pair 
of disk-shaped knots aligned at a distance of 3.49 xlO^'^ cm 
from the central star, along the long axis of the ellipsoid. 
The cylindrical jets widen into cone-shapes at the knot ends 
in order to simulate the effect of material accumulating at 
the knots, as suggested by the HST images (particularly for 
knot K4, as seen in the right panels of Figure 2) . The diam- 
eter of the base of the cones equals that of the disk-shaped 
knots. The centres of the 3.49x10^'' cm diameter circular 
disks representing the knots are aligned with the centres 
of the cylindrical jets (hence they are seen almost edge on). 
The width of the disks is assumed to be 3.88 x 10^^ cm in our 
model, although only a fraction of these is ionised, as is clear 
from the right panels of Figure 2. The H number density in 
the jets and knots is taken to be homogeneous and equal to 
1250 cm"'' and 1500 cm"'', respectively, consistently with 
the values derived in Paper I. 

Results from the ellipsoidal rim and the spherical outer 
shell are combined into a single R-component to enable us 
to carry out a direct comparison with the slit spectra from 
Paper l, as we show in Table 1. 

The jets (J-component) are included in our simulation 
as the radiation field has to be transferred through them 
before reaching the outer ansae {K-component). However, 
given that the emission detected from this region is very 
faint and that such structures may not be in equilibrium, 
we take our results for the J-component as very uncertain 
and omit them from any further discussion. 

Finally, neither the inner caps nor the tenuous halo were 
included in our model. The former because they do not lie o n 
the same axis as the K-component ijReav fc Athertonll98!f) . 
and are therefore not expected to have a major influence on 
the ionization structure of the outer knots. And the latter 
because it is very faint, and therefore it is not expected 
to contribute significantly to the integrated emission line 
spectrum. 

Figure 1 shows model profiles along the major axis, in 
agreement with our assumed density, as Ne/N_r/ ratio is cor- 
related to the level of ionisation of the gas (see Section 3.4). 



3 RESULTS 

The best fit to the observed spectra was obtained assum- 
ing the parameters given in Table 2, as discussed in Sec- 
tion Table 3 shows the predicted and observed intensi- 
ties of some important coUisionally excited lines, in which 
values are given relative to that of H/3(=100) for each neb- 
ular component (R and K) and integrated over the nebula 
(NEB). 

The model H/3 fluxes are given in the first row of the 
table. The dereddened line intensities quoted in the Obs col- 
umn of Table 3 were obtained from the spectroscopic data 
presented in Table 1 of Pap er I, by using a logarit hmic ex- 
tinction constant cna =0.16 llGoncalves et al.ll2003ll and the 
reddening law of ICardelli et aLfi^SsF " For each line inten- 



sity, in each component, the upper row (of the Obs column) 
shows the values for the North-East (Rl and Kl) regions of 
NGC 7009, while those for the South- Western R2 and K4 
components are given in the lower row (see Figure 1 and 
Table 1 of Paper l). 

The line intensities predicted by our model were con- 
volved with a long-slit profile (assumed to be rectangular) 
aligned along the long axis of the PN. This is a necessary 
correction for an extended object with a complex geometry, 
such as that of NGC 7009, if any meaningful conclusions re- 
garding the ionisation and temperature structures are to be 
gained from the comparison of the model with the observa- 
tions. The dimensions of 1.5" vs. 4' at a distance of 0.86 kpc 
were assumed in order to be consistent with Paper I. For 
each nebular component listed in Table 3 we give the slit 
and no-slit line intensities in adjacent columns. 

The absolute value for the observed H/3 flux of 
NGC 7009 of 3.197 x 10"^° erg cm"^ s"\ quoted in Ta- 
ble 3 , was obtained from the VLA radio recombination line 
flux dCarav et alJll989ll . This flux should be compared to 
the nebula-integrated value (NEB) predicted by our model 
for the no-slit case. A very good agreement (better than 2%) 
is found. 



3.1 Narrow slit effects 

From the comparison of the model slit and no-slit columns 
it appears that some emission lines are more affected than 
others. This is very easy to understand if one considers the 
electron temperature distribution and the physical extension 
of the various regions where each of the relevant ionic species 
are most abundant. 

He IIA4686, for example, is enhanced by ~31% in the 
slit results for the R component and ~29% overall. The op- 
posite behaviour is shown by the He I lines for which the slit 
results show a ~5% depletion for He lA5876 and He lA6678 in 
both the R component and NEB component. Similar beav- 
iours are observed for other species as well, particularly we 
note the enhancement of intensities in the slit column for 
[N ll] , [S ll] and [S lll] lines. We also note that different lines 
from the same ionic species appear to be enhanced/depleted 
by different amounts, like [S ll]A4069, A4076, enhanced by 
~10-12% and A6717,A6731 enhanced by -15-18%. This is 
due to the different sensitivities of the various transitions to 
changes in the electron temperatures. 



3.2 Comparison of the emission lines spectrum 

We show in the last three columns of Table 3 a comparison of 
our model with the observations integrated over the whole 
slit. Whilst a satisfactory agreement is obtained for many 
emission lines, some discrepancies do remain, including the 
case of the [N ll] and the [S ll] lines. These discrepancies 
can be readily understood by noticing the different emission 
lines intensities measured for the North-East and the South- 
West sides of the nebula (Table 1 of Paper l). In our models 
we have assumed the nebula to be symmetric about the x-, 
y- and z-axis and cannot therefore reproduce such asym- 
metries, which will be reflected into the integrated nebular 
spectrum for the species that are affected the most. This 
is not a great concern for us, given that a detailed model 
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Table 3. Model and observed dereddened spectra of R, K and NEB. Line intensities are normalized to H/3=100. 



5 



Line Identification (A) R K NEB 

Model Model Obs. Model Model Obs. Model Model Obs. 
no-slit slit no-slit slit no-slit slit 



H/3(10"i^ erg cm-2 g-i) 3119 



397.1 



7.93 



2.94 



3136 



405.7 3197* 



[On] 3726.0 + 3728.8 5.38 5.88 

[Nom] 3868.7 111. 112. 

**[Neni] 3967.5 34.5 34.8 

[Sii] 4068.6 0.31 0.35 

[Sii] 4076.4 0.10 0.11 

[Oiii] 4363.2 8.54 8.55 

Hen 4685.7 16.8 22.0 

**[Ariv] 4711 3.92 3.87 

[Ariv] 4740.2 5.19 5.24 

H/3 4861.3 100. 100. 

[Oiii] 5006.8 1197 1172 

[Cliii] 5517.7 0.35 0.36 

[Cliii] 5537.9 0.57 0.59 

[Nil] 5754.6 0.10 0.11 

Hei 5875.7 14.1 13.4 

[Siii] 6312.1 1.47 1.55 

[Nil] 6583.4 5.87 6.64 

Hei 6678.1 3.99 3.82 

[Sii] 6716.5 0.46 0.53 

[Sii] 6730.8 0.84 0.99 



15.5 
7.34 
106. 
105. 
51.7 
44.6 
1.50 
2.36 
0.87 
1.16 
7.76 
8.65 
26.0 
23.4 
5.48 
4.96 
5.56 
4.54 
100. 
100. 
1162 
1225 
0.43 
0.43 
0.53 
0.55 
0.14 
0.18 
13.9 
14.1 
1.27 
1.14 
7.09 
5.60 
3.87 
3.95 
0.48 
0.38 
0.84 
0.69 



213. 



131. 



40.5 



4.95 



1.63 



10.2 



0.00 



1.50 



1.27 



100. 



1315 



1.09 



1.04 



3.04 



15.8 



4.04 



167. 



4.48 



21.9 



29.3 



241. 
129. 
40.1 
5.46 
1.80 
9.88 
0.00 
1.40 
1.19 
100. 
1272 
1.08 
1.03 
3.48 
15.9 
3.99 
193. 
4.49 
24.2 
32.4 



204. 
155. 
82.3 
141. 
35.3 
59.8 
7.43 
5.04 
2.36 
1.83 
6.94 
9.64 
1.00 
1.23 
2.19 
2.00 
1.14 
1.49 
100. 
100. 
1238 
1310 
0.00 
0.97 
0.00 
0.90 
6.49 
3.91 
18.8 
15.0 
3.89 
3.25 
355. 
194. 
7.49 
3.18 
36.8 
23.2 
51.0 
28.4 



6.06 



111. 



34.6 



0.33 



0.10 



8.55 



16.7 



3.91 



5.17 



100. 



1198 



0.35 



0.57 



0.11 



14.1 



1.48 



6.39 



3.99 



0.53 



0.94 



8.21 
113. 
34.9 
0.40 
0.13 
8.59 
21.6 
3.84 
5.18 
100. 
1177 
0.37 
0.59 
0.15 
13.5 
1.59 
8.44 
3.83 
0.79 
1.32 



24.2 

108. 

48.9 

1.98 

1.11 

8.15 

15.8 

4.10 

3.92 

100. 

1206 

0.54 

0.64 

0.46 

14.5 

1.68 

27. 

3.97 

2.33 

3.85 



The Model spectra are given with and without considering the narrow slit effect. Upper rows give the observed intensities of the 
North-East (Rl and Kl) part of the nebula, while those for the South- West (R2 and K4) zone are given in the lower rows. See Figure 1 
of Paper I for the sizes of Rl, Kl, R2 and K4. "0.00" as the model predicted HellA4686 emission and as the observed [Cllll] AA5517,5737 
from K, former means that model does not produce any such emission and latter means that such lin es were not detecte d in the spectra 
of the knots. *The absolute value of the observed H/3 flux for the whole nebula was obtained from lGarav et al] ll989ft . **From the 
spectra in Paper I what we actually measured was [Nelll]A3967.5+HeA3970.1 and [Ariv] A4711-|-HelA4713. 



specific to NGC 7009 is not being sought, our goal being 
tlie construction of a model able to explain the apparent 
enhancement of some low-ionisation species in LISs of PNe 
such as NGC 7009 in terms of photoionisation effects only, 
without the need of assuming enhancements in the total el- 
emental abundances of those regions. For this reason, in the 
modelling we aimed at reproducing the spectra of the in- 



dividual regions (R and K) with values falling between or 
being close to one of the observational data for North-East 
or the South- West regions. 

A good agreement is shown in the rim and shell (namely 
R, as justified in Section 2.3.2) with nearly all predicted lines 
falling at intermediate values between the two sides of the 
nebula (upper and lower rows of the Obs column), or being 
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Table 4. Mean temperatures (K) weighted by ionic species. 



Ion 



1 J 1 1 . [ 1 iV.l ill 


J 




III 


IV 




VI 


VII 


H 


10,002 
10,625 
10,106 


10,378 
10,562 
10,380 












He 


9,832 
10,603 
9,914 


9,955 
10,562 
9,968 


12,238 
10,565 
12,238 












9,800 


9,874 


10 089 


10 733 


12 887 


12 813 


10 378 




10,628 


10,598 


10,559 


10,506 


10,494 


10,563 


10,563 




10,142 


9,961 


10,101 


10,731 


12,887 


12,813 


10,380 


N 


9,798 


9,875 


10,098 


10,706 


12,757 


13,158 


13,072 




10,684 


10,625 


10,555 


10,507 


10,490 


10,563 


10,563 




10,402 


10,015 


10,110 


10,704 


12,757 


13,158 


13,070 


O 


9,795 


9,863 


10,073 


12,452 


12,962 


13,183 


13,301 




10,697 


10,621 


10,549 


10,515 


10,563 


10,563 


10,563 




10,529 


9,987 


10,082 


12,452 


12,962 


13,183 


13,301 


Ne 


9,805 


9,920 


10,209 


12,590 


13,026 


13,269 


13,312 




10,627 


10,585 


10,561 


10,528 


10,563 


10,563 


10,563 




10,041 


9,967 


10,215 


12,590 


13,026 


13,269 


13,312 


S 


9,793 


9,845 


10,021 


10,423 


11,341 


13,021 


13,308 




10,635 


10,605 


10,560 


10,510 


10,467 


10,469 


10,563 




10,221 


9,973 


10,040 


10,422 


11,339 


13,021 


13,308 



For each element, from the upper to the lower row, we show values for R, K and NEB, respectively. 



within 10-30% of one of the two values. The only notable 
exception is the [S ll]AA4068.6, 4076.4 auroral doublet, with 
both component being depleted by factors of 4 to 10. The 
same behaviour is not observed in the knots. We also note 
that MOCASSIN predicts slightly enhanced values for the [S 
n]AA6717,6731 nebular doublet in the R-component. Since 
the auroral lines have much higher critical densities than 
the nebular lines, it is possible that they are being produced 
in a region that is not accounted for in our model. It would 
be useful to check the behaviour of the [O ll]A7320 doublet, 
also auroral transitions, but unfortunately these fall outside 
the spectral range explored by Paper I. 

The agreement between model and dereddened intensi- 
ties of the K-component is even better than that of the 7?- 
component. All model predictions fall in between the data 
observed at the two sides of the nebula, or are within 10- 
30% of one of the two values. Our model does not predict 
any He llA4686 emission from the knots. We argue that the, 
however quite low, emission reported at the locations of Kl 
and K4 could be attributed to intervening material in the 
PN halo not included in our model. 

In the following we will discuss in more detail the pre- 
dicted mean temperatures and ionic abundances derived 
from our modelling. 

3.3 Mean temperatures 

Mean temperatures weighted by the ionic abundances are 
given in Table 4. Expressions for the volume integrated frac- 
tional ionic abundances (discussed in Section 3.4) as well as 
those for the mean te mperatures wei ghted by ionic abun- 
dances were defined in lErcolano et a l. ( 2003b) , their eqs. 
(2) and (3). The values listed in Table 4 show that the 



temperature structure of the nebula is well reproduced by 
our model. Figures of Te[0 ill] obtained for regions R, K 
and NEB are within the empirical values determined in Pa- 
per I (see Table 1) for the two sides of the nebula, namely 
10,100 K, 10,000 K, and 10,100K for rim, knots and NEB, 
respectively. Similarly, the values of Te[N ll] of the three 
components agree with the empirical values to better than 
5%, while for Te[S II] the agreement is to within 33% and 
11% for the two knots, due to our model underestimating of 
the [S ii] auroral lines (see Section 3.2). 

We confirm that the electron temperature distribu- 
tion is fairly homogeneous acr oss the vari o us re gions 
of the nebula, see fo r instance, iRubin et alJ ll2002h and 
ISabbadin etaU t004) . The mean Te for the whole neb- 
ula, from the neutral (I) to the highly ionized ions (VII), 
for all species, are respectively 10,200±200K, 10,050±150K, 
10,050±850K, 11,400±1050K, 12,600±700K, 13,100±200K 
and 12,650±1,300K. 

Our model results for Ne and Te compare well with 
the many estimates available from the literature; in partic- 
ular Te [O III] and Te [N ll] for the rim, can be compared to 
the values of HyunE & A Ucr ( 1995a b ) ; Ma th is et ah (199j); 
iKwitter fc Henrvl lll998l) : ILuo et al.l feoOllh : ISabbadin et alJ 
J2004h that cluster around 10,000K for Te[0 iii], while val- 
ues between 9,000K and 12,800 K, peaking aroun d 10,550 K, 
were found for TeJN ll] b y 'Hvuue & AUeil lll995allbl) : 
iKwittcr & Henrv ( 1993l : lGoncalves et al. (.2003). As for the 
outer knots, Te[N ll] an d Te[0 lll] were fo und as being 
8,100K a nd 11,500K bv iBalick et alJ (|l993), 9,300K and 
9.1 00K bvlKwitter fc Henrvl (1199811 arid 11,350 and 10,000K 
bv lGoncalves et all tOO^ . 
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Table 5. Averaged fractional ionic abundances. 



Ion 



Element 


I 




11 


in 


IV 


V 


VI 


VII 




H 


5.40 


-4) 


0.999 
















1.04 


-2) 


0.989 
















6.46 


-4) 


0.999 














He 


7.00 


-4) 


0.814 


0.185 














7.39 


-3) 


0.992 
















7.74 


-4) 


0.817 


0.181 












C 


3.10 


-6) 


1.08(-2) 


0.555 


0.428 


5.45(-3) 










2.42 


-4) 


0.132 


0.834 


3.24{-2) 












5.33 


-6) 


1.21(-2) 


0.560 


0.421 


5.33(-3) 








N 


1.37 


-6) 


5.84(-3) 


0.546 


0.443 


4.12(-3) 


1.00(-5) 








3.73 


-4) 


0.136 


0.823 


4.01{-2) 












4.47 


-6) 


7.14(-3) 


0.551 


0.437 


4.03(-3) 


9.85(-6) 






O 


4.83 


-6) 


9.72(-3) 


0.862 


0.125 


3.13(-3) 


l.ll(-6) 








2.86 


-3) 


0.188 


0.808 














2.81 


-5) 


1.15(-2) 


0.863 


0.122 


3.07(-3) 


1.08(-3) 


1.08(- 


6) 


Ne 


3.56 


-6) 


1.07(-2) 


0.917 


7.11{-2) 


7.20(-4) 










1.56 


-4) 


6.67(-2) 


0.933 














4.98 


-6) 


1.14(-2) 


0.918 


6.95(-2) 


7.04(-4) 








S 


6.84 


-7) 


8.13(-3) 


0.343 


0.546 


9.89(-2) 


2.64(-3) 


9.97(- 


6) 




8.16 


-5) 


0.155 


0.755 


8.81{-2) 


6.71(-4) 










1.42 


-6) 


9.71(-3) 


0.350 


0.540 


9.70(-2) 


2.58(-3) 


9.75(- 


6) 


Ar 


1.52 


-7) 


1.25(-3) 


0.363 


0.613 


1.90(-2) 


2.86(-3) 


2.51(- 


5) 




2.99 


-5) 


2.20(-2) 


0.821 


0.156 












4.05 


-7) 


1.47(-3) 


0.370 


0.606 


1.85(-2) 


2.80(-3) 


2.46(- 


5) 



For each clement the first row is for the R, the second is for K, and the last one is for NEB. 



ble to the ionic abundances derived from the observations 
in Paper I, with an agreement better than 90% for aU of 
the above ions, except the S^^/H of R, whose discrepancy 
reaches 70%. 

In Figure we show the ionic abundance profiles of 
the model, as well as the electron density and temperature 
profiles along the axis which includes the knots. The top- 
left panel in this Figure, reflects the Ne (continuous line) 
and Te (dashed line) variations through the long axis of 
the nebula. It shows that Ne peaks at the innermost re- 
gion of the R-component, having a mean value through the 
component of 4,100 cm~^. The J-component as well as the 
K-component have mean Ne of 1385 cm~'^ and 1640 cm"'', 
respectively, while somewhat smoothed proflles are present 
at the edges. The electron density profile is a result of the 
ionization structure combined with the density distribution 
assumed in Section 2.3.3. As for the Te profiles, although 
varying somewhat within the components, have mean val- 
ues that are in close agreement with the values measured 
in Paper I, see numbers quoted in the previous subsection. 
This will be discussed further in Section 4. 

Figure 1 also highlights the strong dependence of the 
ionisation level on the geometry and density distribution of 
the gas. It is therefore clear that an apparent overabundance 
of N+ in the knots can be produced by providing the correct 
gas opacity to screen this region from the direct stellar pho- 
tons. We should add at this point that an even larger N"*" 
abundance can be obtained by further enhancing the gas 
density at the rim-jet interface, without significantly chang- 
ing the [S ll] density ratio in the R-component. 



3.4 Fractional ionic abundances 

Results for the fractional ionic abundances are shown in Ta- 
ble 5. Hydrogen and helium are fully at least singly-ionized 
in R, K and NEB, and significant fractions of the heavy ele- 
ments are in higher ionization stages in the rim as well as in 
the nebula as a whole (NEB). We also note the lower ioniza- 
tion of the knots. An important issue that should be noted 
here is the N/N^ ratio being higher than the O/O^ ratio 
by a factor of 1.39 in the knots, 1.66 rim and 1.61 in the to- 
tal nebula. This result is at variance with the N£Nj_= O / O 
jKinesburgh fc Barlow! 1 19941: iPerinotto et alJl2004b^ gener- 
ally assumed by the ic/ method, with the consequent errors 
on empirically derived total elemental abundances, like those 
in Paper I. 

We note from Table 5, that only a small fraction ( 0.6%, 
14% and 0.7% for R, K and NEB, respectively) of the total 
nitrogen in the nebula is in the form of N" and N"^. As only 
lines from these ions were observed (see Table 3 of Paper l), 
the nitrogen abundance determination is particularly uncer- 
tain. Sulfur abundances also suffer from similar problems. 

Combining the values in Table 5 with those for the to- 
tal elemental abundances used by our model (Table 2) , frac- 
tional ionic abundances relative to H are readily available. 
The total abundance that we have used for all regions mod- 
eled returns ionic abundances, for R and K respectively, of 
He+ = 0.091 and 0.111; N+ = 1.17x10"*' and 2.72x10"''; 

= 3.88x10""* and 3.64x10""*; Ne+^ = 9.72x10"^ and 
9.89x10"^; and finally S+^ = 3.09x10"® and 6.79x10"^ 
These — which are the most important ions of the optical 
range for determining the total abundances — are compara- 
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Distance from the centre (cm) 



Figure 1. Model profiles along the major axis. From left to right and top to bottom, the electron density and temperature, as well 
the ionic fraction for He, C, N, O and S are shown. The position of the R-component (9.9x10^® up to 1.84xl0^^cm) and that of the 
K-component (3.49x10^'' up to 3.87x lO^'^cm) arc marked in each panel. After the inner cavity the inner wall of the ellipsoidal shell is 
encountered, this has Nj:/ = 4,000 cm~^ and Ne ~ 5,000 cm"''. The value of Ne decreases through the first shaded area that represents 
the ellipsoidal shell, as the level of ionisation of the gas decreases. After the first shaded area we find the jet (N^=l,250 cm"'', Ne ~ 
1,500 cm~^). At the end of the jet we find the knot, where and the N^f = 1,500 cm~^, while Ne reaches about 1,700 cm~^. 



3.5 Visualization of the model results 

NGC 7009 has been observed with the HST WFPC2 with fil- 
ters centred in the [O lll]A5007 and [N n]A6583 emission lines 
(see Section 2.1). These archive images, already published in 
Paper I, are compared to the model predicted emission maps 
in Figure |5| The maps were produced at an inclination of 
84° with respect to the line o f sight, as indicated by t he kine- 
matics of the PN polar axis jReav fc AthertoJl985l) . for [O 
III] (right top panel) and [N ll] (lower right panel) emission 
lines. First of all we call the readers attention to the fact 
that the polar bubble, which appear in the HST [O ill] im- 
age, at fainter intensity levels, as an extension of the shell, 
and also the inner pair of knots, K2 and K3 were not con- 
sidered in our modelling. Thus, excluding the polar bubble 
and the inner knots, the images and our maps are at least 
qualitatively in good agreement (note that a photometric 
comparison is beyond the scope of this paper). 

From the maps we note the higher excitation of the 
equatorial axis of the PN as compared to the polar one, be- 
cause of the higher density of the rim, and following the 
pattern known from previous studies, as stated in the Intro- 
duction. One also clearly sees that the [N ll] map is more 
extended than the [O ill] one, as expected from a nebula 
excited by a central star. Finally, one notes that the knots 
are fainter (as compared to the inner regions) in [O ill] and 
brighter in [N ll], due to the enhanced recombination in the 
knots. All the above suggest that our model can satisfacto- 
rily reproduce the main features in the images of NGC 7009, 
that could in principle be extended to other PNe with a 



pair of outer low-ioni zation knots (see, for more examples, 
iGoncalves et alJl200lD . 



4 ALTERNATIVE MODEL 

As noted in previous sections, the ionisation structure of 
the model is strongly dependent on the input 3D density 
distribution. Because of that, here we also present the re- 
sults obtained with a slightly different geometry, which we 
call the alternative model. The only difference in terms of in- 
put density distribution between our initial model and the 
alternative one is that in the latter the jets start 3.3% fur- 
ther away from the central star than in the previous model. 
Since the H number density in the jet is four times lower 
than the value in the rim at the position of the polar axis, 
our change results into an increase in the line optical depth 
along the path connecting the central star to the outer pair 
of knots. The intensity of the ionizing stellar field reaching 
the knots is therefore decreased in our alternative model, 
hence resulting in a lower ionization level. Other model pa- 
rameters were also adjusted in order to optimise our fit of 
the spectrum and these are summarised in Table 6. More 
importantly, with this alternative model we aimed at repro- 
ducing the more extreme [N ll] emission of the Eastern knot 
of the nebula (Knotl in Paper I and in Figure 2). 

As in the case of our original model discussed in Sec- 
tion 3, we found a good agreement of the alternative model 
with the dereddened fluxes of Paper I, for the rim, the knots 
as well as for the whole nebula (NEB). Again, most of the 
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[OIII] shell 



polar bubble 
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Figure 2. Left panel: HST [O III] and [N ll] images of NGC 7009 on a logarithmic intensity scale; slit position for the spectra discussed 
in Section 3 is indicated by short lines (P. A. = 79°), while labels mark the position of the several structures, where "K", "J" and "R" 
stand for "knot", "jet", and "rim", respectively. Right panel: project emission maps from the model; the top map is [O in] and the 
bottom one is [N ll]; maps (in arbitrary units) are not on the same colour intensity scale. 



Table 6. Input parameters for the alternative model. 



model are not shown here as they are qualitatively identical 
to those of our original model (right panels of Figure 2). 



L* (L©) 


3136 


N/H 


2.8x10"* 


Te// (K) 


82,000 


O/H 


4.5x10"* 


Rin (cm) 


0.0 


Ne/H 


1.0x10"* 


Rout (cm) 


3.88x101'^ 


S/H 


0.9xl0"5 


He/H 


0.12 


Ar/H 


1.5x10"*^ 


C/H 


3.5xl0~* 


Fe/H 


1.25x10"® 



Abundances are given by number, relative to H. 

line intensities are within 10-30% of one of the sides of 
NGC 7009, or between the two values of the Eastern and 
Western sides, as quoted in Table 7, with the exceptions 
of the [S n]AA4068,4076 and [S ii]AA6717,6731 doublets, as 
pointed out in Section 3.2. 

The mean Te of the alternative model also compare 
very nicely with the observed values, namely, Te[0 ill] and 
Te[N ii] of 10,200 K and 9,950 K, 10,600 K and 10,600 K, 
and 10,200 K and 10,100 K, for rim, knots and NEB, re- 
spectively. And the Ne are automatically matched with the 
empirical values, since they are constraining the model as 
input parameters. 

The corresponding ionisation structure is shown in Ta- 
ble 8, from which we obtain the N/N"*" and O/O'*' ratios 
of the different zones in the nebula, as being: 1.64 (7?- 
component), 1.33 (K-component) and 1.56 (NEB), that, as 
above, are in contradiction with the ratios adopted in the 
empirical ic/ scheme. We chose to include Table 8 in the pa- 
per, although qualitatively similar to Table 5, to allow for 
model icjs to be derived at a later time should one wish to 
do so. 

Finally, the projected emission maps obtained from this 



5 CONCLUSIONS 

This work focused on the study of the apparent N overabun- 
dance in the outer knots of NGC 7009 with respect to the 

mai n nebular rim and sh ell. 

[Alexander fc Balicl^ ^1997^ and ICruenwald ~ Vieaaj 
showed that long-slit data may give spurious over- 
abundances of N and other elements in the outer regions 
of model PNe. These authors have identified the different 
charge-exchange reaction rates of N and O as the main 
responsible for the effect in the low-ionization regions of 
the nebulae, therefore affecting the (N"'"/N)/(0'''/0) ratio, 
that as discussed in Section 3.4, is commonly u sed to ob- 
tain t he icf for nitrogen. However, as shown by iMampasol 
i2004f) charge-exchange reactions can, in this case, account 
at most by 20% of the nitrogen overabundance of the knots 
of NGC 7009. 

In this paper we have presented a model that was able 
to reproduce the main spectroscopic characteristics of the 
various spatial regions of NGC 7009 without the need of 
assuming an inhomogeneous set of abundances. We investi- 
gated the importance of taking into account the effects of 
a narrow slit, and our results in Table 3 and 7 show that 
the convolution of the model results with the profile of the 
narrow slit used for the observations presented in Paper I 
results in the [N ll] emission being enhanced respect to H/3 
in all regions of the nebula, with the effect being slightly 
more pronounced in the knots. 

The (N+/N)/(0+/0) predicted by our models are 0.60, 
0.72 and 0.62 (or alternatively, 0.61, 0.75 and 0.64) for the 
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Table 7. Altcrnativo model arid olwcrvcd dcix-ddciicd s])(Klra of R. K and NEB. Lino intensities are nornralized to H/3=100. 



Line Identification (A) R K NEB 

Model Model Obs. Model Model Obs. Model Model Obs. 
no-slit slit no-slit slit no-slit slit 



H/3(10-i3 erg cm-^ s"!) 



3120 



398.5 



7.95 



2.95 



3137 



407.1 3197* 



[On] 3726.0 -|- 3728.8 5.37 5.95 

[Nciii] 3868.7 109. 110. 

**[Ncin] 3967.5 33.8 34.0 

[Sii] 4068.6 0.32 0.36 

[Sil] 4076.4 0.10 0.12 

[Oiii] 4363.2 9.11 9.13 

Hen 4685.7 18.4 24.0 

**[Ariv] 4711 4.28 4.23 

[Ariv] 4740.2 5.67 5.73 

H/3 4861.3 100. 100. 

[Oiii] 5006.8 1234 1208 

[Cliii] 5517.7 0.35 0.36 

[Cliii] 5537.9 0.57 0.59 

[Nil] 5754.6 0.15 0.16 

Hex 5875.7 15.0 14.3 

[Sin] 6312.1 1.50 1.57 

[Nil] 6583.4 8.12 9.32 

He: 6678.1 4.25 4.06 

[Sil] 6716.5 0.45 0.54 

[Sii] 6730.8 0.84 1.00 



15.5 
7.34 
106. 
105. 
51.7 
44.6 
1.50 
2.36 
0.87 
1.16 
7.76 
8.65 
26.0 
23.4 
5.48 
4.96 
5.56 
4.54 
100. 
100. 
1162 
1225 
0.43 
0.43 
0.53 
0.55 
0.14 
0.18 
13.9 
14.1 
1.27 
1.14 
7.09 
5.60 
3.87 
3.95 
0.48 
0.38 
0.84 
0.69 



257. 



123. 



38.1 



5.89 



1.95 



9.86 



0.00 



1.49 



1.26 



100. 



1262 



1.07 



1.03 



5.33 



17.0 



3.97 



296. 



4.81 



26.0 



34.9 



299. 
122. 
37.9 
6.72 
2.22 
9.39 
0.00 
1.38 
1.17 
100. 
1205 
1.05 
1.01 
6.32 
17.0 
3.91 
352. 
4.82 
29.8 
39.9 



204. 
155. 
82.3 
141. 
35.3 
59.8 
7.43 
5.04 
2.36 
1.83 
6.94 
9.64 
1.00 
1.23 
2.19 
2.00 
1.14 
1.49 
100. 
100. 
1238 
1310 
0.00 
0.97 
0.00 
0.90 
6.49 
3.91 
18.8 
15.0 
3.89 
3.25 
355. 
194. 
7.49 
3.18 
36.8 
23.2 
51.0 
28.4 



6.18 



109. 



33.8 



0.33 



0.11 



9.12 



18.3 



4.27 



5.65 



100. 



1234 



0.35 



0.57 



0.16 



15.0 



1.50 



9.03 



4.26 



0.54 



0.96 



8.81 
110. 
34.1 
0.42 
0.14 
9.16 
23.5 
4.19 
5.65 
100. 
1212 
0.37 
0.60 
0.22 
14.4 
1.61 
12.5 
4.08 
0.85 
1.41 



24.2 

108. 

48.9 

1.98 

1.11 

8.15 

15.8 

4.10 

3.92 

100. 

1206 

0.54 

0.64 

0.46 

14.5 

1.68 

27. 

3.97 

2.33 

3.85 



As in Table 2, with the top/bottom rows given intensities of the Eastern/ Western size of the nebula. 



rim, knots and the whole nebula, respectively, all at vari- 
ance with the icf assumption of unity for this ratio. The icjs 
will therefore be underestimated by the empirical scheme, 
both in the case of the R and K components, but more 
so in the former (by a factor of 1.21). Therefore this effect 
may partly be responsible for the apparent inhomogeneous 
N abundance derived from observations. The differences in 
the (N"'"/N)/(0"''/0) in these two components may be due 
to a number of effects, including charge exchange, as stated 
above, and the difference in the ionization potentials of the 
relevant species, which makes the (N+/N)/(0^/0) ratio ex- 
tremely sensitive to the shape of the local radiation field. For 



this reason, a realistic density distribution is essential to the 
modelling of a non-spherical object, if useful information 
is to be extracted from spatially resolved observations. The 
density distribution of the gas modifies the shape of the local 
radiation field via the gas opacities and (as demonstrated in 
Section 3.4 and Section 4) matching the emission line spec- 
trum of a given nebular region relies, among other things, 
on the careful construction of a model capable of providing 
the correct optical depths in various directions. 

Our main conclusions discussed here may also be ex- 
tended to other PNe exhibiting FLIERs, such as NGC 6543 



Nitrogen abundance of the NGC 7009 's outer knots 

Table 8. Averaged fractional ionic abundances for the alternative model. 
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Top to bottom rows give values for R, K and NEB, respectively. 



and NGC 6826, for which simi lar apparent chem ical inho- 
mogeneities have been claimed jBalick et alJll994) . 
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